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Abstract 

A  fuel  cell  system  with  planar  stack  configuration  is  studied  with  a  view  of  portable  applications.  A  mathematical  model  to  estimate  the 
discharge  characteristics  of  the  system  is  presented  in  this  research  to  clarify  the  effect  of  the  number  of  segmented  cells  on  the  output  power 
and  energy  losses  generated  in  the  system.  In  the  isothermal  case  of  this  analysis,  we  revealed  that  output  power  of  a  planar  stack  SOFC 
increases  with  increasing  the  number  of  cell  segments  because  the  current  distribution  in  the  electrolyte  becomes  more  uniform  by  increasing 
the  number  and  this  results  in  a  reduction  of  energy  losses  in  the  system. 
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1.  Introduction 

In  the  conventional  fuel  cells,  in  order  to  produce  a  practical 
output  power,  there  is  need  for  stacking  several  cells  with 
bipolar  separators.  The  use  of  separators  make  the  structure  of 
fuel  cells  very  complicated  and  has  been  a  principal  problem 
in  miniaturization  of  the  fuel  cells.  Recently,  we  proposed 
novel  small  size  SOFC  called  SSSC  [1]  in  which  the  seg¬ 
mented  cells  are  connected  in  series  and  placed  in  a  row  in  the 
flow  direction  of  fuel-oxidant  mixture  without  separators.  We 
expect  that  SSSC  (Fig.  1)  is  easier  to  produce  than  that  of 
conventional  SOFC  stack  because  of  no  separator.  The  stack 
configuration  is  called  “planar  stack”  in  this  research,  and  in 
the  previous  paper  [1],  we  dealt  with  an  analysis  of  the  planar 
stack  SOFC  using  very  simple  mathematical  model  which 
employs  Newman’s  assumption  [2]  that  the  sum  of  the 
activation  overpotential  and  concentration  overpotential  is 
constant  on  either  electrodes  over  the  whole  cell  area.  The 
model  is  very  simple  and  if  a  load  current  is  given,  we  can 
obtain  cell  voltage  without  complicated  iteration  procedure, 
however  the  model  is  not  taking  IR  drop  in  both  the  electrodes 
and  potential  distribution  in  the  electrolyte  into  consideration 
and  corresponding  to  the  hypothetical  fuel  cell  modeling 
when  the  electrode  resistances  are  0.  Discharge  characteristics 
of  SOFC  are  strongly  affected  by  not  only  operating  condi¬ 
tions  such  as  cell  current,  temperature  and  partial  pressures 
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but  also  structural  parameters  such  as  thicknesses  and  lengths 
of  electrode  and  electrolyte.  Structural  parameters  are  closely 
related  to  potential  distribution  in  each  component  material, 
and  the  changes  of  the  distributions  result  in  output  power 
change.  Therefore,  for  the  development  of  the  planar  stack 
fuel  cell  having  optimum  dimensions,  it  is  important  to 
estimate  the  influences  of  structural  parameters  on  the  dis¬ 
charge  characteristics  with  more  practical  calculation  model 
considering  potential  distribution  in  each  cell  component. 

In  this  paper,  we  are  dealing  with  fuel/oxidant  separate 
system,  then  fuel-oxidant  mixture  system  is  not  considered. 
Our  aim  is  not  to  analyze  the  performance  of  one-chamber 
fuel  cell.  The  most  emphasis  on  this  paper  is  to  quantita¬ 
tively  clarify  the  effect  of  the  cell  partition  along  gas  flow 
direction  through  the  analysis  of  the  planar  stack  fuel  cells.  If 
the  fuel-oxidant  mixture,  namely  the  partial  oxidation  of 
hydrocarbon  is  considered,  the  system  will  become  very 
complicated.  The  merit  of  the  cell  partition  discussed  in  this 
paper  is  expected  in  all  types  of  fuel  cells,  that  is  to  say  not 
only  in  one-chamber  fuel  cell.  Therefore,  for  the  simplifica¬ 
tion,  we  discussed  only  for  the  fuel/oxidant  separate  system. 

2.  Theory 

For  the  single  cell  schematically  shown  in  Fig.  2,  the 
potential  distribution  can  be  considered  by  the  equivalent 
circuit  of  Fig.  3.  In  the  SOFC,  because  the  resistivity  of  the 
electrolyte  is  considerably  larger  than  that  of  electrode, 
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Nomenclature 

d 

electrolyte  thickness  (m) 

Eloss,  act 

energy  loss  caused  by  activation 
overpotential  (W) 

£loss,  el 

energy  loss  caused  by  IR  drop  in  the 
electrodes  (W) 

£loss,  ion 

energy  loss  caused  by  IR  drop  in  the 
electrolyte  (W) 

emf 

electromotive  force  (V) 

emf* 

net  electromotive  force  (V) 

F 

Faraday  constant  (C  mol-1) 

AG°(7) 

standard  Gibbs  energy  change  at  TK 
(J  mol-1) 

AH°(T) 

standard  enthalpy  change  at  TK  (J mol-1) 

4 

load  current  (A) 

4 

current  in  the  anode  (A) 

4 

current  in  the  cathode  (A) 

4 

current  in  the  electrolyte  (A  m-1) 

k 

rate  constant  (A  m-2) 

L 

cell  length  (m) 

Mi 

molar  flow  rate  (mol  s-1) 

n 

segment  number 

N 

the  number  of  segment  cells 

p 

total  pressure  (atm) 

Pi 

partial  pressure  (atm) 

R 

gas  constant  (J  mol-1  K-1) 

R, 

electric  resistance  of  anode  (Q  m-1) 

Rc 

electric  resistance  of  cathode  (Q  m-1) 

Re 

electric  resistance  of  electrolyte  (Q  m) 

t 

electrode  thickness  (m) 

T 

temperature  (K) 

Uh2 

hydrogen  utilization  (%) 

Wo) 

cell  voltage  at  given  cell  current  (V) 

V,(/0) 

total  cell  voltage  of  planar  stack  fuel 
cell  (V) 

W 

cell  width  (m) 

X 

position  (m) 

X 

molar  flow  rate  of  oxygen  atom  in  the  fuel 
gas  (mol  s-1) 

Greek  letters 

^th 

theoretical  thermal  efficiency  of  fuel  cell 
(%) 

^con,  s 

energy  conversion  efficiency  of  fuel  cell 
system  (%) 

^con,  e 

energy  conversion  efficiency  from  chemical 
energy  to  electricity  (%) 

n 

activation  overpotential  (V) 

r\lR 

IR  drop  (V) 

P 

resistivity  (Q  m) 

Subscripts 

a 

anode 

ave 

average 

c  cathode 

e  electrolyte 

f  fuel  gas  mixture 

max  maximum 

ox  oxidant  gas  mixture 


current  is  assumed  to  pass  parallel  to  gas  flow  direction  in  the 
electrodes  and  normal  to  in  the  electrolyte.  Nagata  et  al.  have 
reported  some  numerical  analyses  on  conventional  SOFC 
characteristics  using  this  model  [3-5].  Based  on  the  model, 
Mathematica™  program  that  can  calculate  partial  pressure, 
current  and  potential  distributions  and  I-V  (current-voltage) 
and  I-P  (current-output  power)  curves  of  a  planar  stack 
SOFC  was  made  in  this  research. 


Fig.  1.  Schematic  diagram  of  a  single  chamber  SOFC  with  series 
connected  cells,  which  is  abbreviated  to  SSSC. 
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Fig.  2.  Single  cell  configuration. 


Emf*(x)=Emf(x)-{ \ric(x)\+T].Jx) } 
Fig.  3.  Equivalent  circuit  of  s  planar  SOFC. 
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In  order  to  understand  the  effect  of  the  cell  partition  along 
gas  flow  direction  intrinsically,  first  of  all,  discharge  char¬ 
acteristics  of  planar  stack  fuel  cell  should  be  analyzed 
neglecting  heat  transfer  in  the  cell.  In  the  case  of  the 
analysis,  temperature  increase  in  the  cell  are  not  taken  into 
consideration,  in  other  words,  all  heats  from  the  cell  caused 
by  IR  drops,  activation  overpotentials  and  reaction  entropy 
are  released  from  the  fuel  cell  system  without  contribution  to 
temperature  increase  in  the  cell.  The  case  of  the  analysis  is 
referred  to  as  “isothermal  case”. 

2.7.  Current  distribution 


In  order  to  consider  the  effect  of  the  activation  over¬ 
potential  on  the  SOFC  output  power,  the  net  electromotive 
force,  emf*(v)  is  introduced.  The  emf*(v)  is  expressed  as  the 
subtraction  of  the  activation  overpotential  {\f]c(x)\  +  f]a(x)} 
from  the  electromotive  force,  emf(x).  By  applying  the 
Kirchhoff’s  law  to  a  small  area  of  v  axis  in  the  model, 
we  have  following  equations: 


h{x)  = 


d/a(x) 


dx 

d{/e(x)7?e(x)} 

dx 


-/c(x)  Rc(x)  +  7a(x)  /?a(x) 
demf  *  (x) 


(i) 


(2) 


Eq.  (1)  means  that  the  difference  in  current  in  the  anode 
Ia(x)  between  v  +  Ax  and  x  is  equal  to  current  in  the  electro¬ 
lyte  Ie(x).  Eq.  (2)  indicates  that  the  current  distribution  in  the 
electrolyte  is  controlled  by  IR  drop  in  the  cathode  and  anode 
and  net  electromotive  force  distribution.  Moreover,  the  fol¬ 
lowing  relationship  holds  among  7a(v),  Ic(x)  and  70: 

7a(*)+7c(*)=70  (3) 


The  differential  equation  related  to  current  in  the  anode 
Ia(x)  is  derived  from  Eqs.  (l)-(3). 


d 2Ia(x) 
dx2 


dIa(x)dRe(x)  1  Rc(x)  +Ra{x) 

dx  dx  Re(x)  aW  Re(x) 


f  demf*  (x) 
{  dx 


1 

Re(x) 


(4) 


We  can  solve  this  equation  by  giving  initial  net  electro¬ 
motive  force  and  some  boundary  conditions,  then  the  current 
distributions  in  the  anode  Ia(x)  can  be  obtained.  When  the 
load  current  is  70,  7a(L)  =  70  and  7a(0)  =  0  are  the  boundary 
conditions.  In  this  case,  7a(L)  =  7o  is  used  to  find  the  value  of 

(dIa(x)/dx)x=o  which  is  essential  for  solving  differential 
equation  of  second  order  (Eq.  (4)).  First,  Eq.  (4)  is  solved 
using  7a(0)  =  0  and  temporary  value  of  (dIa(x)/dx)x=o.  Next, 
the  value  of  (dIa(x)/dx)x=o  is  changed  so  that  7a(L)  =  70  can 
be  satisfied,  and  when  it  is  satisfied,  we  can  obtain  right 

solution.  Using  the  solution,  current  distribution  in  the 
cathode  Ic(x)  and  electrolyte  7e(v)  can  be  obtained  with 
Eqs.  (1)  and  (3),  respectively. 


2.2.  Mass  balance  in  oxidant  gas  mixture 

In  this  paper,  SOFC  for  H2-02  system  is  considered. 
Therefore,  only  oxygen  molecule  is  the  reactive  species  in 
the  oxidant  gas  mixture.  The  cathode  reaction  is  expressed  by: 

|02  +  2e“  — >  02~  (5) 

Hence,  oxygen  atom  transfers  across  the  electrolyte  as  the 
form  of  oxide  ion  from  cathode  side  to  anode  side.  This 
transfer  is  equivalent  to  the  current  flow  in  the  electrolyte 
from  anode  to  cathode.  Then,  oxygen  atom  in  the  oxidant  gas 
mixture  is  consumed  at  any  position  of  v  by  the  quantity 
corresponding  to  Ie(x).  Moreover,  according  to  Eq.  (1), 
current  in  the  anode  at  v  is  equal  to  the  integration  of  current 
in  the  electrolyte  from  0  to  v.  Therefore,  Ia(x)  corresponds  to 
the  amount  of  oxygen  moved  across  the  electrolyte  from  0  to 
v.  Consequently,  the  oxygen  molar  flow  rate  in  the  oxidant 
gas  Mq2  ox  (x)  gradually  decreases  in  the  gas  flow  direction. 
This  behavior  can  be  expressed  by: 

7  (x) 

Mo2,oxW=Mo2,ox(0)-^  (6) 

where  MO2,Ox(0)  is  oxygen  molar  flow  rate  at  oxidant  gas 
entrance. 

2.3.  Mass  balance  in  fuel  gas  mixture 

Hydrogen,  water  vapor  and  nitrogen  gases  are  supplied 
into  anode  side,  therefore  in  the  fuel  gas,  reactive  species  are 
hydrogen,  water  vapor  and  oxygen.  The  anode  reaction  is 
expressed  by: 

H2  +  02 4“  ->  H20  +  2e“  (7) 

Molar  flow  rate  of  oxygen  atom  in  the  fuel  gas  X(x)  is 
expressed  as: 

X(x)=MH2  o,f(0)+7^-  (8) 

where  Mn2o,f(0)  is  molar  flow  rate  of  water  vapor  at  fuel 
inlet. 

Thus,  in  the  fuel  gas  mixture,  we  can  obtain  concentration 
distributions  of  reactive  species  Mq2,  f(x),  MH2,  f(x)  and 
47n2o,fW  from  the  following  simultaneous  equations  [3]: 

7  (x) 

Mn2o,f(x)  +  2Mo2,f  (x)  =  MH2o,f(0)  + 

47n2o,fW  +Mn2jf(x)  =  Mu2  o,f(0)  +MH2,f(0) 

47H2o,f(X)  =  K 

47n2,f(v)  Po2,fW°-5 * 

Eqs.  (9)  and  (10)  indicate  mass  balances  for  the  oxygen  atom 
and  the  hydrogen  atom,  and  Eq.  (11)  indicates  mass  action 
law,  where  Pq2^(x)  is  the  partial  pressure  of  oxygen  and  K 
the  equilibrium  constant  for  the  combustion  of  hydrogen  at 
operating  temperature. 


(9) 

(10) 

(ID 
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2.4.  Potential  distribution 


2.5.  Energy  losses 


Partial  pressure  distribution  of  component  i,  Pfx)  is  given 
by: 


Pito 


Mi 

E  Mi 


P 


(12) 


where  M  is  molar  flow  rate  and  P  the  total  pressure.  In  the 
oxidant  gas  mixture,  components  are  oxygen  and  nitrogen, 
in  the  fuel  gas  mixture,  hydrogen,  water  vapor,  oxygen  and 
nitrogen. 

Distribution  of  electromotive  force  can  be  calculated 
by: 


where  indices  ox  and  f  means  oxidant  and  fuel,  respectively. 

When  the  dissociation  and  adsorption  of  the  reactive 
species  is  rate  determining  step,  Butler-Volmer  equations 
for  the  current  in  the  electrolyte  are  expressed  as  following 
equations  [6]: 

for  the  cathode  reaction 


Ie(x)  =  WkcPlf0X(x) 


(14) 


for  the  anode  reaction 


exp  'fkf ,,a  "  exp  {  kf ,,a  (x) } 

(15) 

where  kc  and  ka  are  the  rate  constants  of  the  cathodic  and 
anodic  reactions,  respectively. 

As  the  partial  pressure  distributions  and  current  distribu¬ 
tion  in  the  electrolyte  have  been  obtained,  distributions  of 
the  activation  overpotentials  r]c(x)  and  r]a(x)  can  be  calcu¬ 
lated  using  Eqs.  (14)  and  (15).  Moreover,  by  subtracting  the 
anode  and  cathode  activation  overpotentials  from  electro¬ 
motive  force  as  expressed  in  Eq.  (16),  we  can  obtain  net 
electromotive  force  distribution  emf*(v). 

emf*  (x)  =  emf(x)  -  {\t]c(x)\  +  f/a(x)}  (16) 

Eq.  (4)  can  be  solved  again  using  this  new  emf*(x),  and 
then  the  new  distributions  of  current,  concentration, 
partial  pressure,  emf  and  activation  overpotential  can 
be  obtained.  These  calculations  are  iterated  until  conver¬ 
gence  is  achieved.  After  the  convergence  is  achieved, 
if  the  IR  drops  along  current  path  is  subtracted  from 
emf*(v),  cell  voltage  at  given  cell  current  V(/0)  can  be 
obtained  [3]. 

V(/0)  =  emf* (x)  -  {t]IR  a(x)  +  r]IR  t{x)  +  r]IR  c(x)}  (17) 

where  r]/R,  aW?  Oir,  eW  and  rjIR t  c(x)  are  IR  drops  in  the 
anode,  electrolyte  and  cathode,  respectively. 


h{x)  =  -WKP\if  (x) 


In  the  equivalent  circuit  described  in  Fig.  3,  the  energy 
losses  for  the  activation  overpotentials  £iOSs,  act*  f°r  the  IR 
drop  in  the  electrodes  Eloss ?  el  and  for  the  IR  drop  in  the 


electrolyte  £loss?  ion  are  expressed  by  following  equations: 

£loss,  act  — 

[  h{x){\r]c(x)\  +  rja(x)}dx 

Jo 

(18) 

£)os,el  = 

/W-(x,<k+ 

/o  a  4 w  Jo  cV  ;  tcW 

(19) 

^loss,ion  = 

fL  d 

Jo  llix)ps{x)—6x 

(20) 

2.6.  Energy  conversion  efficiency 


Theoretical  thermal  efficiency  of  fuel  cell  eth  is  given  by: 


A  G°(T) 
A//°(298  K) 


x  100 


(21) 


where  AG°(7)  is  standard  Gibbs  energy  change  of  hydrogen 
combustion  at  T  K  and  A#°(298  K)  is  lower-heating  value 
(LHV)  of  hydrogen  combustion. 

Energy  conversion  efficiency  of  fuel  cell  system,  as  can  be 
expressed  by: 


W/o) 

Mh2(0)  A//°(298  K) 


x  100 


(22) 


where  Vt(/0)  is  total  cell  voltage  of  planar  stack  fuel  cell.  On 
the  other  hand,  energy  conversion  efficiency  from  chemical 
energy  to  electricity  se  can  be  expressed  by  following 
equation: 


£e 


W(!o) 

Gh2Mh2(0)  Atf°(298K) 


x  100 


(23) 


2.7.  Calculation  of  the  output  power  for  planar 
stack  fuel  cells 

As  shown  in  Fig.  4,  in  the  planar  stack  fuel  cell,  single  cell 
is  divided  into  several  segments  and  they  are  placed  in  the 
direction  of  the  gas  flow.  The  number  of  segments  is  denoted 
by  N  and  each  cell  is  written  by  segment  number  n.  In  the 
five  segments  cell,  n  =  1  denotes  first  segment  cell  namely 
gas  inlet  side  cell,  on  the  other  hand  n  =  5  denotes  fifth 
segment  cell,  namely,  gas  outlet  side  cell.  For  the  simplifica¬ 
tion,  interconnection  and  the  effects  of  the  mixing  of  the 
cathode  and  anode  side  gases  in  the  gap  between  the 
neighbor  cells  are  neglected.  First,  I-V  curve  for  the  first 
segment  is  calculated  by  above  mentioned  procedure  and  the 
values  of  partial  pressures  at  gas  exit  of  the  first  segment  are 
obtained  for  a  given  load  current  /0.  Then,  I-V  curve  for  the 
second  segment  cell  can  be  calculated  by  the  same  procedure 
using  these  partial  pressure  values.  The  total  cell  voltage  and 
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Fig.  4.  Configuration  of  the  planar  stack  fuel  cell. 


electric  power  for  a  planar  stack  SOFC  are  obtained  as  the 
sum  of  the  values  of  all  segments. 

The  total  cell  length  and  cell  width,  L  and  W,  are  fixed  at  3 
and  1  cm  because  these  dimensions  are  suitable  for  portable 
electric  devices,  and  N  is  varied  from  1  to  5.  The  structural 
parameters  and  operating  conditions  are  summarized  in 
Table  1 .  The  conditions  of  Table  1  are  based  on  the  experi¬ 
mental  conditions  in  the  measurement  of  the  discharge 
characteristics  of  single  chamber  SOFC  reported  by  Hibino 
et  al.  [7]. 


3.  Results  and  discussion 


£/h2  is  introduced  as  an  operating  condition  in  place  of  cell 
current  /0,  because  this  parameter  is  more  useful  than  70  to 
compare  the  discharge  properties  of  single  cell  and  planar 
stack  cell.  Uu2  is  expressed  as: 


Un2  = 


Mo 

2FMH2,f(0) 


(24) 


The  relationship  between  Uu2  and  70  at  present  conditions 
is  shown  in  Fig.  5.  Cell  current  in  the  N  segments  cell 
becomes  1  IN  of  single  cell.  However,  current  density 


Table  1 

Values  used  in  the  calculation 


Fig.  5.  Relationship  between  hydrogen  utilization  and  cell  current  applied 
in  planar  stack  SOFC.  Calculation  conditions  are  equal  to  Table  1. 


expressed  with  I0/(LW/N)  (dotted  line)  is  determined  by 
UK2  independently  of  the  number  of  segmented  cells. 

3.1.  Output  power 

Fig.  6  shows  the  relationship  between  hydrogen  utiliza¬ 
tion  and  output  power  of  planar  stack  SOFC  calculated  using 
the  model  of  Fig.  3.  It  is  obvious  that  single  cell  with  length 
of  3  cm  generates  only  50  mW  at  the  most.  However,  output 
power  of  SOFC  remarkably  increases  by  dividing  a  single 
cell  into  segment  cells  in  the  gas  flow  direction.  In  the  five 
segments  cell,  more  than  200  mW  can  be  generated  at 
Un2  =  80%  (105  mA  cm-2).  This  value  is  considered  to 
be  useful  for  portable  electric  devices.  These  results  are 
summarized  in  Table  2.  In  the  following  section,  the  reason 
for  the  output  power  increase  by  cell  partition  will  be 
discussed  from  the  viewpoint  of  energy  losses. 

3.2.  Energy  losses 

Fig.  7  shows  the  energy  losses  for  the  activation  over¬ 
potentials  £ioss,  act?  IR  drop  in  the  electrolyte  £iOSs,  ion  and  IR 


Structural  parameters 

Cell  width,  W  (cm)  1 

Cathode  thickness,  tc  (cm)  0.0017 

Anode  thickness,  ta  (cm)  0.0017 

Electrolyte  thickness,  d  (cm)  0.05 

Operating  conditions 

Temperature,  T  (K)  1223 

Oxygen  partial  pressure  at  oxidant  inlet,  Pq2  (0)  (atm)  0.165 

Hydrogen  partial  pressure  at  fuel  inlet,  Pu2  (0)  (atm)  0.082 

Water  vapor  partial  pressure  at  fuel  inlet,  Fh2o(0)  (atm)  0.058 
Flow  rate  of  inlet  gases,  /  (cm-3  s-1)  2.5 

Other  parameters 

Resistivity  of  LSM,  pc  (Q  cm)  0.0083 

Resistivity  of  Ni,  pa  (Q  cm)  0.001 

Resistivity  of  YSZ,  pe  (Q  cm)  13.8 

Rate  constant  for  the  cathodic  reaction,  kc  (A  cm-2)  0.218 

Rate  constant  for  the  anodic  reaction,  ka  (A  cm-2)  0.427 


Fig.  6.  Output  power  of  planar  stack  SOFC  for  various  segment  numbers. 
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Table  2 

Dependence  of  discharge  characteristics  of  planar  stack  SOFC  on  the  number  of  segment  cells  for  isothermal  case 


Output  power  (mW) 

£loss,  el 

(mW) 

£loss, 

ion  (mW) 

£loss,  act  (mW) 

£S  (%) 

UH2  =  50% 

Uh2  =  80% 

UKl  = 

50%  Un2  =  80% 

Un2  -- 

=  50%  Un2  =  80% 

Uu2  = 

50%  Uh2  =  80% 

Un2  =  50% 

{/h2  =  80% 

Single  SOFC 

N=  1 

23 

102 

25 

20 

4.9 

Planar  stack  SOFC 

N=  2 

112 

114 

34 

87 

13 

35 

12 

30 

23 

23 

N=  3 

133 

168 

18 

46 

11 

27 

9.5 

25 

27 

34 

N  =  4 

142 

190 

11 

29 

9.6 

25 

8.7 

23 

29 

38 

N=  5 

146 

201 

7.7 

20 

9.3 

24 

8.4 

22 

30 

41 

drop  in  the  electrodes  7qoss ?  el,  respectively.  As  seen  in  this 
figure,  since  current  must  be  applied  to  thin  electrode  layers, 
terrible  energy  loss  takes  place  in  the  electrodes  of  single 
cell  due  to  excessively  long  length  in  the  current  flow 
direction.  In  the  present  conditions,  7qoss?  ei  is  considerably 
larger  than  E\ oss?  act  and  7qoss?  ion  when  the  number  of  seg¬ 
ment  cells  are  less  than  3,  however,  in  the  five  segments  cell 
7qoss?  ei  becomes  smaller  than  the  other  energy  losses.  Fig.  8 
shows  the  relationship  between  cell  current  70  and  7qoss?  ei  f°r 


The  number  of  segment  cells  N 


Fig.  7.  Dependence  of  energy  losses  in  the  planar  stack  SOFC  on  the 
number  of  segment  cells. 


Fig.  8.  Relationship  between  cell  current  applied  in  planar  stack  SOFC 
and  energy  loss  caused  by  1R  drop  in  the  both  electrodes. 


various  segment  numbers  N.  It  is  clear  that  7qoss?  ei  in  planar 
stack  SOFC  is  approximately  determined  by  load  current  70. 
Then,  because  cell  current  of  planar  stack  SOFC  becomes  1/ 
N  of  single  cell  at  a  same  £7h2  as  shown  in  Fig.  5,  7qoss?  ei can 
be  dramatically  reduced  by  increasing  the  number  of  seg¬ 
ment  cells  N.  On  the  other  hand,  as  shown  in  Fig.  7,  7qoss?  act 
and  7qoss ;  ion  can  also  be  reduced  by  increasing  N.  As  can  be 
seen  in  Fig.  9(a),  in  the  single  cell,  current  in  the  electrolyte 
significantly  increases  near  the  gas  inlet  side  for  the  follow¬ 
ing  two  reasons.  First  is  the  difference  in  fuel  concentration 
between  gas  inlet  and  outlet.  As  shown  in  Fig.  10(a),  fuel 
concentration  decreases  towards  downstream,  therefore, 
emf(v)  decreases  along  gas  flow  direction.  This  results  in 
current  concentration  at  gas  inlet  side  due  to  highest  driving 
force  for  electrochemical  reaction.  This  effect  is  considered 
at  the  third  term  in  the  right  side  of  Eq.  (2).  Second  is  the 
difference  in  electrode  resistances  between  cathode  and 
anode.  This  effect  corresponds  to  the  first  and  second  terms 


Fig.  9.  Distribution  of  current  in  the  electrolyte  at  Uu2  =  50%:  (a)  for  the 
single  SOFC  and  (b)  for  the  planar  stack  SOFC  (three  segments  cell). 
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(a)  Fuel  decrease 
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ft 
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Fuel  decrease 


Emf(U)  t 
x  =  L 


(b)  Electrode  resistance 


in  the  right  side  of  Eq.  (2).  When  the  cathode  and  anode 
resistances  do  not  agree  with  each  other,  cell  current  will 
pass  in  favor  of  in  the  electrode  that  has  lower  resistance.  At 
present  conditions,  LSM  has  about  10  times  as  large  resis¬ 
tance  as  Ni  has.  Therefore,  current  passes  mainly  in  Ni  as 
shown  in  Fig.  10(b).  This  is  the  main  cause  of  current 
localization  at  gas  inlet  side.  Naturally,  current  concentra¬ 
tion  observed  in  Fig.  9(a)  gives  rise  to  activation  over¬ 
potential  increase  nearby  the  gas  inlet  side  on  both  the 
electrodes.  Concentrations  of  current  and  activation  over¬ 
potential  at  gas  inlet  side  are  the  reasons  why  7qoss t  act  and 
Eioss,  ion  of  the  single  cell  become  large  as  seen  in  Fig.  7. 
However  it  is  clear  that  these  energy  losses  decrease  with 
increasing  N ,  because  as  shown  in  Fig.  9(b)  current  in  the 
electrolyte  becomes  more  uniform  by  dividing  a  single  cell 
into  segments.  Current  concentration  in  the  electrolyte  is 
defined  as  the  ratio  of  maximum  current  7e  max  to  average 
current  7e  ave.  By  the  way,  according  to  Eq.  (2),  if  the 
resistance  parameters  of  cathode  and  anode,  Rc(x )  and 
Ra(x ),  are  set  at  0,  we  can  obtain  current  distribution  in 
the  electrolyte  caused  by  fuel  decrease  in  the  gas  flow 
direction,  therefore  current  concentration  due  to  fuel  dis¬ 
tribution  also  can  be  estimated.  The  relationship  between  the 
number  of  segment  cells  and  current  concentration  is  shown 
in  Fig.  11.  As  for  planar  stack  SOFC,  7e  max  is  the  average 
value  of  that  in  each  segment  cell.  As  observed  in  Fig.  1 1,  the 
current  concentration  decreases  with  increasing  N  with  the 
relation  of  7e?  max/7e  ave  =  6.0A-0'9,  and  it  is  clear  that 
current  concentration  in  the  electrolyte  of  the  single  cell 
is  caused  by  mainly  the  difference  in  electric  resistance 
between  cathode  and  anode.  For  the  single  cell,  80%  of  the 
current  concentration  is  due  to  the  difference  in  the  electric 
resistances.  In  contrast,  for  the  five  segments  cell,  80%  of 


T=1223  K 

LXW=3  cmX  1  cm 

O  Total 

A  Caused  by  fuel  decrease 

Black  Un2  =  80  % 

White  Uh2  =  50  % 

\y  =  6.0jt<>9 

y=  Oy^ 

.  . . . 

The  number  of  segment  cells  N 


Fig.  11.  The  effect  of  the  number  of  segment  cells  on  current 
concentration  in  the  electrolyte. 


that  is  due  to  the  fuel  decrease  in  the  gas  flow  direction. 
Dependence  of  7qOSs,  act  and  7qoss?  ion  on  the  current  concen¬ 
tration  is  shown  in  Fig.  12.  As  can  be  seen  in  Fig.  12,  by 
increasing  N ,  at  a  same  hydrogen  utilization  or  current 
density,  7qoss?  act  and  7qoss ?  ion  can  be  reduced  depending 
on  the  decrease  of  the  current  concentration  in  the  electro¬ 
lyte,  and  the  change  of  the  energy  loss  value  Alog  7qoss  is 
determined  by  the  change  of  the  current  concentration 
A(7e  max/7e?  ave).  These  reductions  of  energy  losses  by  cell 
partition  correspond  to  the  output  power  increase  shown  in 
Fig.  6.  The  dependence  of  energy  losses  on  the  number  of 
segment  cells  is  also  summarized  in  Table  2. 

3.3.  Energy  conversion  efficiency 

Dependence  of  energy  conversion  efficiency  on  the  num¬ 
ber  of  segment  cells  is  shown  in  Fig.  13.  The  £s  is  very  small 
when  N  is  equal  to  1  because  of  large  IR  drop  in  the 
electrodes  and  its  maximum  value  is  about  10%  when 
£7h2  is  30  %.  However,  if  N  and  Uu2  are  increased  up  to 
5  and  80%,  respectively,  system  efficiency  increases  up  to 


Fig.  12.  Dependence  of  energy  losses  caused  by  activation  overpotentials 
Eioss,  act  and  IR  drop  in  the  electrolyte  £ioss,  ion  on  the  current 
concentration. 
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Hydrogen  utilization  /  % 

Fig.  13.  Energy  conversion  efficiency  of  planar  stack  SOFC  for  various 
segment  numbers. 


around  40%.  On  the  other  hand,  £e  increases  up  to  near  the 
theoretical  value  by  decreasing  Uu2  down  to  near  0%. 
Especially  in  the  five  segments  cell,  ee  is  70-50%  in  the 
Un2  range  from  10  to  80%.  Fig.  13  shows  that  cell  partition 
which  leads  to  remarkable  reduction  of  energy  losses  as 
discussed  in  Section  3.2  is  an  effective  way  of  fuel  cell 
efficiency  increasing. 

4.  Conclusions 

The  configuration  of  planar  stack  fuel  cell  is  most  pro¬ 
mising  design  for  portable  electric  devices,  however  so  far 
its  discharge  characteristics  have  not  been  analyzed  quanti¬ 
tatively.  Therefore,  prior  to  fabrication  of  this  fuel  cell 
system,  it  is  very  important  to  clarify  its  discharge  char¬ 
acteristics  using  mathematical  model.  In  order  to  analyze  the 
properties  of  planar  stack  fuel  cell  in  SOFC  field,  Mathe- 
matica™  program  to  calculate  I-V  curve  of  planar  stack 
SOFC  taking  electric  resistances  of  component  materials 
and  activation  overpotentials  into  consideration  has  been 
made. 


Single  cell  of  small  size  SOFC  can  not  generate  sufficient 
electric  power  to  run  portable  electric  devices.  However, 
output  power  can  be  increased  at  a  certain  hydrogen  utiliza¬ 
tion  by  dividing  a  single  cell  into  segment  cells  along  gas 
flow  direction,  in  other  words,  employing  planar  stack 
configuration  because  energy  losses  for  the  IR  drops  in  both 
the  electrodes,  IR  drop  in  the  electrolyte  and  the  activation 
overpotentials  can  be  considerably  reduced  by  cell  partition. 
In  particular,  energy  losses  for  IR  drop  in  the  electrolyte  and 
the  activation  overpotentials  are  closely  related  to  current 
concentration  in  the  electrolyte,  and  the  reduction  of  current 
concentration  by  cell  partition  is  one  of  the  reasons  for  high 
electric  power  of  planar  stack  SOFC. 

Only  hydrogen/oxygen  separate  system  was  considered  in 
this  paper  for  simplification.  But  the  dealt  system  is  not 
suitable  for  the  electric  power  source  of  portable  electronics 
in  the  points  of  operating  temperature  and  considered  fuel. 
The  one-chamber  SOFC  developed  by  Hibino  et  al.  which 
can  be  operated  at  low  temperatures  below  773  K  [8]  is  most 
promising  system  in  the  view  of  portable  applications.  In  the 
SOFC  field,  the  planar  stack  configuration  will  become  more 
effective  for  practical  applications,  especially  when  it  is 
incorporated  into  the  one-chamber  fuel  cells  because  of  no 
need  of  separators  (Fig.  1). 
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